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The  damage  threshold  for  brittle  materials  Impacted  by  water  drops  has  been 
studied  by  using  a combination  of  analytic  and  numerical  approaches  to  obtain 
solutions  for  the  Impact  stress  fields  and  dynamic  stress  Intensity  factors. 

It  has  been  shown  that  expressions  for  the  damge  threshold,  which  consists  of 
the  activation  of  pre-existent  surface  microcracks,  can  be  derived  fi^  these 
solutions.  The  expressions  Identify  three  types  of  behavior,  depending  upon 
the  adjacence  of  the  pre-existent  cracks  to  the  Impact  center,  vis-a-yls  tteir 
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size,  and  describe  the  respective  roles  of  the  three  Important  target  parameters 
the  fracture  toughness,  the  elastic  wave  velocity  and  the  pre-existent  flow 
size.  An  absolute  damage  threshold  has  also  been  defined.  The  conditions  at 
this  threshold  are  demonstrated  to  be  reasonably  consistent  with  available 
threshold  measurements. 
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SUMMARY 


The  damage  threshold  for  brittle  materials  Impacted  by  water  drops 
has  been  studied  by  using  a combination  of  analytic  and  numerical  approaches 
to  obtain  solutions  for  the  Impact  stress  fields  and  dynamic  stress  Intensity 
factors.  It  has  been  shown  that  expressions  for  the  damage  threshold,  which 
consists  of  the  activation  of  pre-existent  surface  microcracks,  can  be 
derived  from  these  solutions.  The  expressions  Identify  three  types  of  be- 
havior, depending  upon  the  adjacence  of  the  pre-existent  cracks  to  the 
Impact  center,  vis-a-vis  their  size,  and  describe  the  respective  roles  of 
the  three  Important  target  parameters,  the  fracture  toughness,  the  elastic 
wave  velocity  and  the  pre-existent  flaw  size.  An  absolute  damage  threshold 
has  also  been  defined.  The  conditions  at  this  threshold  are  demonstrated 
to  be  reasonably  consistent  with  available  threshold  measurements. 


I . INTRODUCTION 
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Th«  damage  threshold  In  brittle  materials  (particularly  ceramic 

materials)  Impacted  by  water  drops  consists  of  the  propagation  of  pre- 

7 ? 3 

existent  surface  cracks  to  form  well-defined  circumferential  cracks  ’’ 

(Fig.  1).  The  cracks  are  activated  by  a tensile  pulse  associated  with  the 

Rayleigh  wave  that  emanates  from  the  Impact  site.  The  Initiation  of 

2 

damage  typically  occurs  In  the  subsonic  regime,  and  has  been  shovm  to 
depend  primarily  on  three  target  properties:  the  fracture  toughness,  the 
elastic  wave  velocity  (tlnodul us/densi ty ] and  the  sl2e  distribution  of  pre- 
existent surface  cracks;  while  the  Important  projectile  properties  are, 
evidently.  Its  size  and  velocity.  Damage  parameters  have  been  derived 
that  Indicate  the  trend  toward  an  Increasing  threshold  velocity  as  the 
target  toughness  or  wave  velocity  Increase.  However,  these  parameters 
have  not  yet  been  sufficiently  quantified  to  permit  the  effective  prediction 
of  threshold  velocities.  In  addition,  the  role  of  the  Initial  flaw  size 
distribution  remains  ambiguous. 

A quantitative  definition  of  the  damage  threshold  will  ultimately 
require  an  experimental  Investigation.  This  study  will  Inevitably  entail 
a comprehensive  series  of  Impact  tests:  searching  for  the  onset  of  damage. 
Before  embarking  on  such  a program  some  theoretical  Insights  regarding 
the  expected  trends  would  be  of  considerable  benefit.  This  would  permit 
the  efficient  choice  of  test  conditions  and  of  materials  that  probe  the 
most  critical  trends.  The  Intent  of  the  present  paper  Is  to  develop  the 
requisite  quantitative  approach  to  the  damage  threshold  problem.  The 
concept  Involves  the  use  of  selected  numerical  calculations  as  specific 
Inputs  to  crack  activation  and  growth  analyses:  ultimately  leading  to 
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to  the  development  of  expressions  for  the  threshold  velocity. 

The  numerical  Inputs  required  for  analysis  of  the  damage  threshold 

are  the  dynamic  stresses  produced  by  the  Impact,  and  the  time  dependent 

2 

stress  Intensity  factors  at  pre-existent  cracts  . Appropriate  solutions 
are  obtained  by  analyzing  the  response  of  ZnSe  Impacted  by  water  drops. 
Normalized  forms,  suitable  'for  more  general  use,  are  then  described.  There* 
after,  an  analytic  approach  is  used  to  obtain  expressions  for  the  damage 
threshold,  by  Invoking  approximate  analytic  functions  for  the  stress  fields 
derived  from  the  numerical  computlons.  The  stress  Intensity  factor  func- 
tions thus  derived  are  compared  with  the  specific  numerical  solutions  to 
substantiate  the  results.  The  final  expressions  for  the  threshold  are 
obtained  In  terms  of  the  critical  velocity  for  damage  Initiation,  as  a 
function  of  the  Important  physical  properties  of  the  target  and  water  drop. 

2.  NUMERICAL  IMPACT  SOLUTIONS 
2.1  The  Numerical  Method 

The  calculations  are  conducted  using  a finite-difference  computer 
code,  WAVE-L.  This  code  Is  a two-dimensional  (axl-synmetrlc  or  planar), 

4 

Lagranglan,  explicit  code  based  on  the  HEMP  scheme  . It  Integrates  the 
governing  partial  differential  equations  of  motion  for  arbitrary  dynamic 

problems  concerned  with  solids  or  fluids.  It  has  been  applied  to  a number 

5 6 

of  Impact  Investigations,  under  both  subsonic  and  hypersonic  conditions. 

The  method  employed  In  the  present  Investigation  Is  to  use  a 
pressure  loading  function  previously  obtained^  by  examining  the  pressures 
developed  during  Impact  of  a rigid  surface  by  a water  drop  (F1g.  2).  This 
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function  1s  pertinent  whenever  the  surface  deflection  Induced  by  Impact  Is 
negligible  compared  with  the  drop  radius.  It  Is  thus  a reasonable  ap» 
proximatlon  for  high  modulus  target  materials  at  low  (e.g..  subsonic) 

[ ' Impact  velocities.  This  Is  the  condition  of  present  Interest.  The  pre- 

sure  function  p(r,t)  specifies  both  the  spatial  and  temporal  features  of 
the  Impact  pressure. 

The  pressure  function  derived  for  a water  drop  with  a radius,  rp, 
of  1 mm  Impacting  a rigid  surface  at  a velocity  of  220  ms~^  Is  used  to 
calculate  the  stresses  that  develop  In  a ZnSe  target  material.  The  results 

Q 

P can  be  normalized  using  the  dimensionless  parameters  described  by  Blowers 

I (Eqn.  16). 

i 

I 2.2  Impact  Stresses 

I 

Some  preliminary  calculations  were  performed  to  examine  the  magni- 
tude of  the  near  surface  stress  gradient,  by  conducting  stress  calculations 
I at  a specific  location  (the  position  on  the  surface  where  the  tensile  stress 

Is  a maximum)  as  a function  of  cell  size.  The  results,  plotted  In  Fig.  3, 
demonstrate  that  the  stress  Increases  rapidly  as  the  surface  Is  approched: 

I O 

a result  anticipated  by  Blowers’  analysis  . This  feature  will  become  an 

Important  aspect  of  the  damage  threshold  model.  Because  of  the  rapid 

I 

spatial  variation  In  stress,  all  subsequent  studies  are  conducted  using 
I the  finest  grid. 

^ The  dynamic  stress  field  results  of  primary  Interest  for  damage 

! analyses  are  the  temporal  trends  In  the  stress  (particularly  the  radial 

I tensile  stresses  responsible  for  circumferential  cracks)  at  various  spatial 

I locations,  and  the  spatial  stress  variations  normal  and  parallel  to  the 
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surface  at  times  approaching  the  development  of  the  peak  tensile  stress 
(at  that  location]. 

Initially,  In  order  to  provide  a useful  perspective  of  the  stress 
fields,  loci  of  the  peak  radial  tensile  stress  (I.e.,  the  stress  regardless 
of  time)  have  been  computed  as  a function  of  the  distance  from  the  Impact 
center,  at  several  locations  beneath  the  surface  (Fig.  4).  Near  the  surface, 
the  peak  tension  develops  at  a radial  distance,  r = r ■ Q^Zr^*  But,  the 
tension  decays  appreciably  with  depth  and  the  peak  value  tends  to  displace 
toward  larger  values  of  r. 

The  stress  gradients  are  of  more  specific  Interest  for  damage  cal- 

2 

culatlons  . The  stresses  computed  at  several  radial  locations,  at  the 
Instant  where  the  surface  tension  attains  Its  maximum  value,  are  plotted 
In  Fig.  5.  The  gradient  diminishes  as  r Increases,  as  anticipated  from 
the  previous  results.  The  stress  distribution  can  be  approximately  ex- 
pressed by  an  exponential  relation  of  the  form; 


o-vo  (r)  exp  [ - (z/z.)  ] 


(1) 


where  a Is  the  peak  stress  at  the  surface  and  Is  a scale  dimension  that 

Increases  as  r Increases;  the  values  of  a and  z^  are  tabulated  In  Table  I. 

g 

Note  that  the  exponential  spatial  dependence  Is  typical  of  Rayleigh  waves  . 

The  second  stress  field  characteristic  of  Interest  Is  the  shape 
of  the  stress  pulse^^.  The  tensile  pulse  computations  are  summarized  In 
Fig.  6.  The  primary  tensile  pulse  Is  located  at  the  Rayleigh  wave  front, 
as  noted  previously.  The  pulse  decreases  In  amplitude  and  tends  to  broaden 
as  the  wave  propagates  outward.  The  broadening  Is  possibly  attributed  to 
varying  contributions  from  the  longitudinal  wave^^:  as  suggested  by  the 
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tendency  for  a low  amplitude  tension  to  be  retained  to  larger  times  at 
the  more  remote  locations.  The  primary  tensile  pulses  are  relatively 
.symmetric  at  all  locations.  The  peak  stress  amplitude  tends  to  decrease 
as  the  distance  from  the  surface  Increases,  but  occurs  at  approximately 
the  same  time  (given  the  radial  location)  for  all  depths  at  which  an 
appreciable  tension  Is  experienced.  The  pulses  also  tend  to  broaden 
slightly  as  the  depth  Increases. 


2.3  Stress  Intensity  Factors 

Stress  Intensity  factor  calculations  are  performed  by  Inserting 
cracked  elements^*^  (I.e.,  elements  that  do  not  support  tensile  stresses 
or  In-plane  shear  stresses]  at  specific  locations  and  computing  the  stress 
gradient  In  the  elements  ahead  of  the  crack.  The  stress  Intensity  factor  K* 
can  then  be  derived  by  fitting  the  crack  tip  stress  field,  taken  along  the 
crack  plane,  to  the  relation^ ^ 


o 


K 

L2itx]^ 


+ 


o“ 


(2) 


where  x Is  the  distance  from  the  crack  tip  and  a®  Is  the  constant  first 
term  of  a Taylor  series  expansion  of  the  non-singular  stress  field  compo- 
nents near  the  crack  tlp^^.  The  procedure  Is  facilitated  by  rearranging 
Eqn.  (2)  Into  the  form, 

o(2irx)**  ■ o®(2itx)'*  + K (3) 

Then,  by  plotting  o(2itx)**  vs  x**,  K can  be  obtained  directly  from  the  Inter- 
cept at  X**  • 0.  The  method  Is  Illustrated  In  Fig.  7:  used  to  obtain  the 
peak  value  of  the  stress  Intensity  factor  K for  each  of  four  crack  lengths, 
n 

All  K values  refer  to  mode  I,  or  Ki,  values;  because  cracks  In  brittle 
materials  under  tension  are  usually  observed  to  follow  mode  I trajectories. 
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at  r/r  ■ 1.5.  It  Is  evident  from  the  linearity  at  small  x that  the  method 

Is  quite  appropriate  for  approximate  K computations;  except  for  the 

smallest  crack  length  (a/rp  ■ 0.0125),  where  the  crack  has  a dimension  of 

just  one  element.  The  method  is  now  used  to  obtain  the  time  dependence 

of  the  stress  intensity  factor  (excluding  the  smallest  crack  length)  for 

cracks  at  r/r  ■ 1.5.  The  results  are  plotted  on  Fig.  8,  The  first  feature 

to  note  is  that  the  stress  Intensity  factors  are  larger  for  the  smaller 

cracks.  The  other  important  characteristic  is  that  K reaches  Its  peak 

value  at  larger  times  for  the  larger  cracks,  even  though  the  peak  tension 

along  the  prospective  crack  plane  occurs  at  approximately  the  same  Instant 

(Fig.  6).  Both  of  these  features  are  explained  in  the  following  section. 

The  former  behavior  Is  related  to  the  stress  gradient  acting  on  the  crack, 

and  the  latter  to  the  time  taken  to  develop  a stress  intensification  at 

the  crack  tip. 

Finally,  the  stresses  computed  at  locations  remote  from  the  crack 
are  used  to  establish  the  magnitude  of  the  near  surface  stress  relief 
created  by  the  crack.  The  results,  plotted  on  Fig.  9,  show  that  the  stress 
near  the  surface  is  appreciably  reduced;  thereby  providing  a rationale  for 
the  spacing  of  the  circumferential  cracks  observed  after  water  drop  impact 
Note,  however,  that  the  stresses  are  enhanced  at  greater  depths:  suggesting 
the  potential  for  crack  instabilities. 

3.  CRACK  FORMATION  THRESHOLD 

The  incidence  of  damage  in  brittle  materials  appears  to  be  adequately 

described  by  the  requirement  that  the  stress  intensity  factor  K at  pre- 

12 

existent  cracks  reaches  the  critical  value  K'  . The  damage  threshold  is 

V 
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analyzed  by  examining  the  development  of  the  stress  intensity  factor  at 
surface  cracks  and  determining  the  conditions  under  which  K reaches 

3.1  The  Stress  Intensity  Factor 


Analyses  of  the  stress  intensity  factor  at  finite  stationary  cracks 

13 

subjected  to  a normal  tensile  stress  pulse  have  indicated  that  K varies 
with  time,  t;  initially  as  t**  and  then  as  a damped  sinusoid  about  the  quasi- 
static stress  intensity  factor  (Fig.  10).  The  initial  motion  of  a 
crack  subjected  to  such  a rapidly  oscillating  stress  intensity  factor  pro- 
bably occurs  by  means  of  small  advance  and  arrest  Increments,  prior  to 
continuous  extension.  The  details  of  this  behavior  are  Ignored  in  the 
present  analysis;  Instead,  It  is  considered  that  the  initial  motion  of  the 
crack  will  be  continuous,  and  described  by  a monotonically  increasing  stress 
intensity  factor  (Fig,  10).  The  function  chosen  to  represent  the  time  and 
crack  length,  a,  dependence  of  K (plotted  in  Fig.  10)  is; 


(K/Kj)^  . [1  + ^2/4)  (a/tc)]-''  (4) 


where  c is  the  wave  velocity  (the  Rayleigh  velocity  for  surface  cracks). 
For  time  dependent  pulse  profiles,  a solution  for  K can  be  obtained  by 
sunning  K due  tJ  each  prior  time  increment,^® 


KCt„) 


(5) 


The  quasi -static  stress  intensity  factor  can  be  obtained  from  the  stress 
o{z,r)  associated  with  the  Rayleigh  wave  by  applying  the  Green's  func- 
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2(L)^  C qCz.r)  [1  •►Ffz/andz 

^ J (a^  - 
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(6) 


where 


FC2/a)  * (1-z/a)  [0.295(z/a)^  + 0.77(z/a)^  - 0.99(z/a)® 


+ 0.51 (z/a)®] 


The  stress  distribution  otz,r)  obtained  from  the  numerical  results 
(Eqn.  CT))  can  be  Inserted  into  Eqn.  (6)  to  give, 

Kj  - 23  ^ Ka/z^j)  (7) 


where  Ua/z^)  Is  plotted  In  Fig.  11:  a convenient  analytic  approximation 
for  Ka/z^),  Is  found  to  be; 


I » Tr/2  '''i+TaTz^  C8) 

so  that 

Kj  = SCira/O+a/Zo)]’*  i aS  (9) 

This  relation  for  can  be  used  in  conjunction  with  Eqn.  (5)  to  obtain 
K(tQ)  if  the  time  variation  of  o is  prescribed.  The  pulse  profiles  pre- 
dicted by  the  numerical  calculations  (Fig.  6)  can  be  approximated  by 
triangular  pulses,  characterized  in  terms  of  stressing  rate  and  an 
unstressing  rate  o^.  Substituting  these  ratios  into  Eqs.  (5)  and  (9) 
yields  the  stress  intensity  factor 
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The  ntinerical  results  CFIg.  6)  indicate  that  a^/Oy **  1 . Applying  this 

ratio  of  stress  rates,  K/^^ci  derived  from  Eqn.  OD  is  plotted  in  Fig.  12, 

as  a function  of  t7t„,  for  several  8.  It  is  noted  that  the  peak  stress 
0 m 

A 

intensity  factor  K develops  after  the  maximum  stress  has  been  reached, 
and  that  the  time  differential  increases  as  the  crack  length  Increases. 

These  characteristics  are  in  accord  with  the  numerical  results.  An  approxi- 

A 

mate  expression  for  K obtained  from  Fig.  12  is; 

= 5 la I*  (T2) 

(1+a/Zg)  £1  + (ir^/4)  (aojj/co^)]  J 


Note  that  the  quantity,  aoj^/cOjj,,  is  just  the  ratio  of  the  time  taken  for  a 
stress  wave  to  traverse  the  length  of  the  crack,  t^(=a/c)  to  the  loading 
time  tj^.  The  numerical  results  (Section  2)  have  been  obtained  for  the  con- 
ditions a > z . t 2t(>,  whereuDon  Eqn.  (12)  reduces  to  the  aoproximate 

O’  C X 

form: 


in 


K - 2(S„)' 
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The  trend  for  R to  decrease  as  the  crack  length  Increases  was  one  of  the 
principal  characteristics  of  the  numerical  results.  A more  detailed  com- 
parison based  on  Eqn.  (12)  is  shown  in  Fig.  13.  The  quantitative  correspon- 
der.ce  is  quite  reasonable  considering  the  approximate  nature  of  both  the 
calculations  and  the  analytic  functions.  We  now  procede  to  develop 

A 

expressions  for  the  threshold,  by  equating  K to  the  critical  stress  Intensity 
factor,  K^. 

3.2  The  Critical  Crack  Extension  Stress 

The  value  of  the  peak  stress  at  the  threshold,  o^,  can  be  found 

A 

t;'  equating  K In  Eqn.  (12)  to  K^,  to  obtain; 


K. 

o • — = cos(*/3) 

/¥T 


(14a) 


where 


♦ ■ cos* 


7.50j^  a' 


KcC/1  + a/Zg 


Two  limit  solutions  are  of  interest.  For  long  relative  loading  times 
(co^  » aojj). 


(14b) 


whereas,  for  short  relative  loading  times  (ca  « aa^^); 


0 ,r i 

C - ■ 


(14c) 
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The  effective  loading  time  will  tend  to  vary  with  both  the  Impact  con- 
dition and  the  distance  from  the  Impact  site.  At  sites  remote  from  the 
Impact  center  the  loading  times  tend  to  be  relatively  long  and  Eqn.  (14b) 
will  tend  to  pertain:  Indicating  that  the  critical  stress  will  either 
be  Independent  of  crack  length  (a  » z^)  or  decrease  according  to  the 
familiar  Inverse  square  root  dependence  (a  < z^).  However,  close  to  the 
Impact  site,  where  the  stress  levels  are  larger,  the  loading  times  are 
relatively  small  and  Eqn.  (14c)  will  tend  to  be  more  pertinent:  sug- 
gesting a regime  In  which  the  critical  stress  Increases  as  the  crack 
length  Increases.  The  Influence  of  the  crack  length  Is  thus  sensitively 
dependent  on  the  location,  v1s-a*v1s  the  Impact  site,  and  on  the  Impact 
condition. 


For  a material  subjected  to  coninuous  Impact  by  water  drops, 
each  crack  In  the  material  will  be  exposed  to  wide  range  of  loading  times 
and  stress  gradients.  A detailed  analysis  of  the  threshold  will  thus 
entail  a determination  of  the  minimum  possible  o^  (for  the  specified 
range  of  Impact  conditions)  and  Its  dependence  on  the  crack  length. 

However,  It  may  also  be  of  merit  to  recognize  the  existence  of  an  absolute 
minimum  In  the  threshold  (as  observed  In  other  threshold  problems)^^'^^. 

An  absolute  threshold  pertains  whenever  a maximun  exists  In  the  variation 
of  the  stress  Intensity  factor  with  crack  length,  as  In  the  present  problem 

A 

(Fig.  13).  The  crack  length  a^  at  the  maximum  stress  Intensity  factor 
Is  obtained  directly  by  differentiation  of  Eqn.  (12)  as. 
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Insertion  of  this  crack  length  Into  Eqn.  (12)  yields  an  expression  for 
the  minimum  possible  stress  that  can  be  tolerated.  Independent  of  crack 
length.  The  resultant  expression  Is  umrieldy,  and  can  only  be  used  to 
full  advantage  when  functional  relations  for  z^Cr]  and  Oj^Cr]  become  avail- 
able. Further  discussion  of  this  threshold  Is  deferred  to  the  following 
section. 


3.3  The  Threshold  Velocity 


The  threshold  can  be  expressed  In  terms  of  the  critical  target 
and  projectile  parameters  by  Invoking  the  dimensionless  material  Independent 

Q 

quantities  T,  Z,  P,  utilized  by  Blowers  In  his  analytic  stress  analysis, 
viz. 


T 


2c!t  2(cz). 

— 5-  , Z - i . n 

(rvlp  (rv)p 


(16) 


where  the  subscripts  p and  t refer  to  properties  of  the  projectile  and 
target  respectively.  Inserting  these  quantities  Into  Eqn.  04)  we  obtain 
a threshold  Impact  velocity  Vp.  It  Is  found  that  three  conditions  prevail, 
determined  by  the  relative  values  of  the  quantity  ; ■ rpVp/ac  and  the 
two  material  Independent  quantities  T^^  (the  dimensionless  time  taken  to 
attain  the  peak  stress  at  the  location  of  Interest  vis-a-vis  the  Impact 
center)  and  Zq  ( a dimensionless  measure  of  the  sub-surface  stress  gra- 
dient at  the  same  location).  The  three  results  are;  for  c larger  than 
both  2/T|_  and  2/1^, 
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for  2/Tj^<  C < 2/Zq, 


1/3  9 1/3 


f i 


vj  - (2/irZ^n^)  ' (K^)  (rcV)p 


1/3 


(17b) 
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I 


1. 

I 

i 


I 


and  for  c<2/T|^,  2/Z^, 

_ A.  9 1/4  O ■5  1/4  _l/9 

vj  ' (w/Z^Tl  n (KgC^a)^  (rcp)p'/2  (17c) 

A 

Where  the  quantity  n Is  the  dimensionless  peak  tensile  stress  at  the  flaw  i 

location.  The  dimensionless  quantities,  Tj^,  Z^  and  12,  can  be  obtained  by  | 

numerical  computation.  Some  typical  values  derived  from  the  results  pre-  | 

sented  In  Section  2 are  sunaarized  In  Table  II,  expressed  as  a function  of 
the  dimensionless  distance  R(  ■ 2(cr)^/(rv)p)  from  the  Impact  center.  Highly 
approximate  functional  relations  for  these  quantities.  In  the  range  R>7 
(the  range  where  fractures  are  usually  observed  ),  derived  from  the  present 
computations  are; 


n > 2.2  (1  - R/22) 


(18) 


Zp  ' 0.08R 

More  precise  functional  relations  can  be  evidently  developed  by  a compre 
hensive  series  of  numerical  computations. 
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These  functional  expressions  for  the  dimensionless  Impact  stress 
field  quantities  can  be  used  to  obtain  a relation  for  the  absolute  damage 
threshold.  This  Is  achieved  by  Inserting  the  crack  size  a^  at  the  threshold 

A A 

Eqn.  (15)  Into  the  general  expression  for  K (Eqn.  12),  equating  K to 

and  then  expressing  the  result  In  terms  of  the  dimensionless  variables 

(Eqn.  16).  The  final  result  Is; 

.2 


. C.3 


(1?) 


■ 

where  X Is  a material  Independent  constant.  The  approximate  magnitude  of 

A 

X can  by  found  by  Inserting  the  quantities  Tj^,  and  from  Eqn.  (18) 
Into  Eqn.  (19),  differentiating  with  respect  to  R and  setting  to  zero  for 
the  mlnimun.  This  procedure  yields  the  result; 


R*  - 10,  X - 1.4 

where  R*  Is  the  value  of  R that  corresponds  to  the  threshold.  It  Is  In- 
teresting to  note  that  Eqn.  (19)  predicts  an  absolute  damage  threshold  for 
Impact  with  1 irm  radius  water  drops  of  128  ms'^  for  ZnS  and  504  ms“^  for 
hot  pressed  SI^N^.  The  former  Is  within  the  range  of  experimental  obser- 
vation^^, providing  encouragement  for  further  development  of  the  present 
approach. 

The  physical  location  of  the  crack  at  the  threshold  condition  R* 

Is  also  of  Interest.  This  location  Is  very  close  In  each  Instance  to  the 
location  of  the  maximum  crack  length  for  ZnS  Impacted  by  250  urn  radius 
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Nylon  spheres  (Table  III).  This  comparison  suggests  that  the  specified 
location  is  indeed  a preferred  site  for  crack  extension,  and  provides 
further  credence  to  the  approach. 

4.0  DISCUSSION  AND  CONCLUSIONS 

The  damage  threshold  associated  with  the  impact  of  a brittle 
surface  by  water  drops  has  indicated  that  the  threshold  velocity  v^  has 
three  forms,  depending  on  the  crack  location  v1s>a>vis  the  impact  center, 
the  crack  size  and  the  impact  condition.  The  conditions  under  which  each 
form  of  the  threshold  should  be  expected  to  pertain  depend  upon  the  rela* 
tive  magnitude  of  the  ratio,  5,  of  the  product  of  the  projectile  radius  "nd 
projectile  velocity  to  the  product  of  the  crack  length  and  the  elastic  wave 
velocity  in  the  target.  The  specific  magnitudes  of  c required  to  invoke 
each  regime  relate  to  the  magnitudes  of  the  material  Independent  impact 
stress  field  variables  T^  and  Z^.  Taking  the  approximate  values  for  these 
variables  obtained  from  the  present,  limited  set  of  numerical  computations 
(Eqn.  18),  the  following  conclusions  can  be  established. 

A 'quasi'Static'  regime  exists  whenever  both  rp>0.07  a(c^/Vp) 
and  r^  > 12.5a,  where  r^  is  the  distance  of  the  crack  from  the  Impact 
center.  This  behaviour  thus  pertains  when  the  cracks  are  small  with  respect 
to  both  the  projectile  radius  and  the  distance  from  the  impact  center 
and  when  the  projectile  velocity  is  large  enough  with  respect  to  the 
elastic  wave  velocity  that  relatively  large  absolute  values  of  the  pulse 
duration  can  be  established*.  The  crack  then  behaves  as  if  it  were  in 

*This  velocity  requirement  derives  from  Blowers' result  for  normalized  time 
(Eqn.  16).  This  normalization  Ignores  some  of  the  complexities  of  the 
contact  process  and  may  not,  upon  more  detailed  examination, be  a stringent 
requirement. 
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a uniform,  quas1>stat1c  field.  Under  these  conditions,  the  threshold 
velocity  Increases  as  the  target  toughness  Increases  and  the  flaw  size 
decreases  (I.e. , as  K^/ZT;  typical  of  quasi-static,  uniform  stress  fields), 
and  decreases  as  the  projectile  density  or  wave  velocity  Increase. 

At  the  opposite  extreme,  a 'dynamic'  regime  exists  whenever  both 
rp  < 0.7aCc^/Vp)  and  r^  < 12.5a.  This  behavior  prevails  when  the  pulse 
duration  Is  small  enough  with  respect  to  the  stress  wave  transit  time  across 
the  crack  that  the  dynamic  response  prevails.  This  requires  relatively 
large  cracks,  close  to  the  Impact  center.  For  such  conditions,  the  threshold 
velocity  Increases  as  the  toughness  and  wave  speed  In  the  target  Increase, 
and  also  Increases  as  the  flaw  size  Increases.  Additionally,  the  threshold 
decreases  as  the  projectile  density,  radius  and  wave  speed  Increase. 

For  Intermediate  conditions,  where  only  one  of  the  above  Inequallo 
ties  are  satisfied,  the  threshold  Is  Independent  of  the  crack  size,  but 
otherwise  depends  on  the  same  valables  as  the  threshold  In  the  dynamic 
regime. 


The  general  trend  In  the  threshold  with  crack  size  thus  exhibits 
Increasing,  constant  and  decreasing  characteristics.  This  is  typical  be- 
havior for  cracks  subjected  to  stress  gradients^ and  leads  to  the  concept 
of  an  absolute  threshold.  The  absolute  threshold  pertains  whenever  a broad 
size  range  of  pre-existent  micro-cracks  exist,  such  that  a crack  always 
exists  that  yields  the  maximum  possible  stress  intensity  factor  for  each 
Impact  situation.  A large  number  of  impacts  are  needed  to  observe  this 
threshold.  The  critical  velocity  at  the  absolute  threshold  depends  on  two 
target  properties,  and  c^,  according  to  It  also  depends  on 
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three  projectile  properties,  rp,  Pp  and  Cp,  In  the  form,  rp  ' pp  ^ Cp  . 

A proportionality  constant  relates  the  threshold  velocity  to  these  target 
and  projectile  properties.  An  approximate  value  for  this  proportionality 
of  1.4  has  been  derived  from  the  present  limited  set  of  numerical  computa- 
tions. 

The  predicted  dependence  of  the  threshold  on  the  target  properties, 
has  been  previously  proposed  as  a damage  parameter^,  and  shown 
to  be  reasonably  consistent  with  some  sparse  experimental  data  on  crack 
activation  by  Nylon  sphere  Impacts.  However,  Interpretation  difficulties 
concerned  with  asigning  the  appropriate  value  for  K were  encountered, 
especially  with  coarse  grained  materials  In  which  the  pre-existent  micro- 
crack lengths  are  on  the  order  of  the  grain  size.  The  Importance  of 

and  c^  and  the  Implications  for  fabricating  materials  with  an  enhanced 

2 

damage  tolerance  have  been  discussed  at  length  In  a previous  paper  . Here 

It  Is  merely  re-emphasized  that  a fine-scale  microstructure  Is  essential, 

and  that  toughening  mechanisms  that  do  not  Invoke  a large  process  zone  must 

be  sought;  an  Important  toughening  mechanism  In  this  category  Is  the  mar- 

19  20 

tensitic  toughening  process  * . Surface  compressloi!  stresses  should  also 

prove  beneficial:  by  reducing  the  Intensity  of  the  Impact  stresses  In  the 
near  surface. 

Absolute  predictions  of  the  threshold  velocity  Indicate  values  of 
128  ms"^  for  ZnS  and  504  ms"^  for  SI^N^.  The  former  seems  reasonable  on 
the  basis  of  the  available  damage  data  for  ZnS,  and  encourages  further 
development  of  this  approach.  The  considerable  superiority  of  SI^N^  Is  also 
apparent,  primarily  by  virtue  of  Its  relatively  high  toughness. 


18 


Hodnid  Intinwaonil 
SC5023.nFR 


The  solutions  for  the  threshold  velocity  now  appear  to  be  suf- 
ficiently quantitative  that  critical  experiments  can  be  designed  to  Investi- 
gate the  threshold.  In  particular,  projectile  Impact  conditions  close 
to  the  predicted  threshold  can  be  selected  to  critically  examine  effects 
of  mlcrostructure.  The  microstructural  effects  can  then  be  Interpreted 
through  their  Influence  on  the  near-surface  fracture  toughness  and  elastic 
wave  velocity:  coupled  with  Influences  of  the  pre-existent  crack  size  distri- 
bution when  the  absolute  threshold  condition  Is  exceeded. 
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TABLE  1 

STRESS  GRADIENT  PARAMETERS 

r/r  oCMPa)  z^/r^ 

1 500  0.015 

1.5  340  0.025 

2 230  0.030 

TABLE  II 

typical  values  of  the  material  independent 

IMPACT  PARAMCTERS 


R 

A 

Q 

^0 

7.5 

1.43 

0.56 

11.2 

0.97 

0.93 

14.9 

0.66 

1.11 
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THE  LOCATION  OF  THE  MAXIMUM  CRACK  LENGTH 
FOR  ZnS  IMPACTED  BY  NYLON  SPHERES^ ^ 


Velocity  (ms”b 

r^(um) 

R*  • 2c^r^/rpVp 

406 

100 

8 

633 

200 

10 

1000 

400 

13 

Calculated  R*'x<10 
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The  pressures  developed  during  laipact  of  a water  drop  onto 
a rigid  surface. 


The  stress  distribution  beneath  the  surface  at  the  location 
of  peak  near  surface  tension*  obtained  for  three  grid 
dimensions. 


Loci  of  the  peak  tenslona,  obtained  for  several  planes 
beneath  the  surface. 
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The  stresses  coaputed  at  a radial  location  r/r  ■ 2 for  a 
crack  positioned  at  r/f  ■ 1.5. 


Fig.  13:  A comparison  of  the  variations  of  the  peak  stress  intensity 

factor  K with  crack  length  a predicted  by  the  analysis  with 
the  values  determined  by  numerical  computation. 


